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The effect of an external magnetic field in QCD phase diagram, namely, 
in the the location of the critical end point (CEP) is investigated. Using the 
2-1-1 flavor Nambu-Jona-Lasinio model with Polyakov loop, it is shown that 
when an external magnetic field is applied its effect on the CEP depends on 
the strength of the coupling. If the coupling depends on the magnetic field, 
allowing for inverse magnetic catalysis, the CEP moves to lower chemical 
potentials eventually disappearing, and the chiral restoration phase transi¬ 
tion is always of first order. 

PACS numbers: 24.10.Jv, ll.lO.-z, 25.75.Nq 

1. Inverse magnetic catalysis in the PNJL model 

The influence of strong external magnetic fields on the structure of the 
QCD phase diagram is a very important held of research due to its conse¬ 
quences on several physical phenomena: the measurements in heavy ion col¬ 
lisions at very high energies, the behavior of the hrst stages of the Universe 
and the understanding of compact astrophysical objects like magnetars. 

The inclusion of a magnetic held in the Lagrangian density of the Nambu- 
Jona-Lasinio (NJL) model and of the Polyakov-Nambu-Jona-Lasinio (PNJL) 
model gives rise to the Magnetic Catalysis (MC) effect, i.e., the enhance¬ 
ment of the quark condensate due to the magnetic held inisii, but fails 
to account for the Inverse Magnetic Catalysis (IMC) found in LQCD calcu¬ 
lations [H El E] where the suppression of the quark condensate takes place 
due to the strong screening effect of the gluon interactions. In order to 
overcome this discrepancy, it was proposed, by using the SU(2) NJL model 
[7] and the SU(3) NJL/PNJL models [8], that the model coupling, Gg, can 
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Fig. 1. (Left panel) The renormalized critical temperatures of the chiral transition 
(T^{eB = 0) = 178 MeV) as a function of eB in the NJL model with a magnetic 
field dependent coupling GgieB) and a constant coupling and the lattice results 
[3]. (Right panel) The chiral {T^) and deconfinement (T)f) transitions temperatures 
as a function of eB in the PNJL, using Gs{eB) given by Eq. ((T|). 


be seen as proportional to the running coupling, and consequently, a 
decreasing function of the magnetic field strength allowing to include the 
impact of as{eB) in both models. Indeed, the strong screening effect of the 
gluon interactions in the region of low momenta weakens the interaction 
which is reflected into a decrease of the scalar coupling with the intensity 
of the magnetic field [9]. 

Since there is no LQCD data available for as(eB), by using the NJL 
model we can fit Gs{eB) in order to reproduce the pseudocritical chiral 
transition temperatures, Tc{eB), obtained in LQCD calculations [3]. The 
resulting fit function that reproduces the Tc{eB) is 


a(C) 


/l + aC^ + 6C^\ 

\i + cc^ + dcy 


( 1 ) 


with a = 0.0108805, b = -1.0133 x 10-^ c = 0.02228, and d = 1.84558 x 
10“^ and where C = We also have used Aqcd = 300 MeV. 

In the NJL model, the renormalized pseudocritical chiral transition tem¬ 
peratures, Tc/Tc{eB = 0), are plotted in left panel of Fig. [T]as a function 
of eB: with the magnetic field dependent coupling Gs{eB) (green line), 
given by Eq. ©; with LQCD results (red dots); and the usual constant 
coupling Gs = G^g (black dashed dot line), that shows magnetic catalyzes 
with increasing Tc /Tc{eB = 0) for all range of magnetic fields. 

Now, using Gs{eB) given in Eq. ([1]), we calculate the chiral and decon- 
hnement transitions temperatures as a function of eB in the PNJL model. 
The results are shown in the right panel of Fig. [1) due to the existing cou¬ 
pling between the Polyakov loop field and quarks within the PNJL model, 
the Gs{eB) does not only affect the chiral transition but also the decon- 
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Fig. 2. Location of the CEP on temperature vs baryonic chemical potential p,B 
(left) and temperature vs baryonic density ps (right) diagrams, for Case 1. The 
baryonic density pB is in units of nuclear saturation density, po = 0.16 fm“^. 


finement transition. Consequently, both temperatures transitions decrease 
with increasing magnetic hied strength. 

2. The influence of the inverse magnetic catalysis in the location 

of the critical end point 

The nature of the phase transition and the existence of the critical end 
point (CEP) are open issues for theoretical studies about the QCD phase 
diagram [10] . From the experimental point of view the existence/location 
of the CEP is also a very timely topic. This renders important to know the 
conditions that can change the position of the CEP in the phase diagram, 
namely the presence of strong magnetic helds. 

In the following, we will study two scenarios for the effect of a static 
external magnetic field on the location of the CEP when symmetric matter 
{flu = fJ^d = fJ^s) is considered: 

Case I — where we take the usual Gg = and no IMC effects are included; 
Case II — where we will use Gs{eB) given by Eq. ([1]) which will allow us to 
consider the IMC effects on the QCD phase diagram. 

The results for Case I are plotted in the left panel of Fig. Oand reproduce 
qualitatively the results previously obtained within the NJL model in m- 
as the intensity of the magnetic field increases, the transition temperature 
increases and the baryonic chemical potential decreases until the critical 
value eB ~ 0.4 GeV^. For stronger magnetic fields both T and fiB increase. 
In the right panel of Fig. [2] the CEP is given in a T versus baryonic density 
plot. It is seen that when eB increases from 0 to I GeV^ the baryonic 
density at the CEP increases from 2pQ to 14/90 m- 

With respect to Case II the results for the CEP are presented in Fig. 
El red points. We clearly observe a different behavior when compared with 
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Fig. 3. Location of the CEP on temperature vs baryonic chemical potential /rs 
(left) and temperature vs baryonic density pB (right) diagrams, for both cases. 
The baryonic density pB is in units of nuclear saturation density, po = 0.16 fm“^. 
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Fig. 4. Masses of the quarks as function of pB for the respective for both 

Cases. 


Case I (black points): at B = 0 both CEP’s coincide but, already for 
small values of B, the CEP is moved to lower temperatures and chemical 
potentials. Nevertheless, until eB ~ 0.3 GeV^ the pattern is similar for 
both Cases. However, for stronger magnetic fields the position of the CEP 
in Case II oscillates between T ~ 169 and T ~ 177 MeV while the chemical 
potential takes increasingly smaller values: a completely different behavior 
when compared with Case I, where both values of T and pB for the CEP 
increase. 

The reason of this behavior lies in fact that the restoration of chiral sym¬ 
metry is stressed by the decreasing of the coupling Gs{eB). The increasing 
of the magnetic filed is not sufficient to counteract this effect as can be seen 
if Fig. m where we plot the quarks masses (M^-black line; M^-red line; Mg- 
blue line) as function of pB for the respective temperature where the CEP 
occurs at eH = 0.1 and eB = 0.5 GeV^. At eH = 0.1 GeV^ (left 

panel) Gg is barely affected by the magnetic held, the values of the quark 
masses are very close to each other for both cases, and the CEP occurs at 
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smaller temperatures and at close, but smaller, chemical potentials. When 
eB = 0.5 GeV^, the quark masses in Case I are increased with respect to 
the B = D case (due to MC effect), being the restoration of chiral symmetry 
more difficult to achieve. However, when Gs{eB), Case II, the masses of 
the quarks are already smaller than the B = 0 case (due to IMC effect) 
leading to an faster restoration of chiral symmetry at small temperatures 
and chemical potentials. Eventually, with the increase of B, the CEP would 
disappear in the temperature axis and the transition to the chiral restored 
phase is always of first order. 
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